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Abstract. Following the tentative identification of the blue luminescence in the Red Rectangle bv lViihetalJ l2005). we com- 
pute absolute fluxes for the vibrational IR emission and phosphorescence bands of three small polycyclic aromatic hydrocar- 
bons. The calculated IR spectra are compared with available ISO observations. A subset of the emission bands are predicted to 
be observable using presently available facilities, and can be used for an immediate, independent, discriminating test on their 
alleged presence in this well-known astronomical object. 
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1. Introduction 



Recentrv. lyiih et alJ J2004 detected a blue luminescence (BL) 
in the w ell-known bipolar protoplanetary Red Rectangle (RR) 
nebula dCohenlll975t ICohen et aljEool . Shortly thereafter, 
IViih et alJ J2005I) . after a thorough analysis of the extinction 
properties and of the correlation of the BL with the observed 
near-IR emission at 3.3 /mi, concluded that they are both most 
likely due to small, neutral polycyclic aromatic hydrocarbons 
(PAHs), composed of three to four aromatic units. Another pos- 
sible explanation of the BL tentatively identifi es it as fluore s- 
cence by ultrasmall silicon nanoparticles dNavfeh et alJ2005T) . 

The hypothesis of the ubiquitous presence of free 
gas-phase PAHs in the interstellar medium (ISM) o riginated 
about 20 years ago dLe ger & Pugedfl984l lAlla mandola etal] 
Il985h . Due to their spectral properties and their high photo- 
stability, these molecules were suggested as the most natu- 
ral interpretation for the so-called "Aromatic Infrared Bands" 
(AIBs), a set of emission bands observed near 3.3, 6.2, 7.7, 
8.6, 11.3 and 12.7 Ami, in many dusty environments excit ed by 
UV photons JLeger et alJll989HAllamandola et alJll989l) . The 
AIBs are the spectral fingerprint of t he excitation of vibration s 
in aromatic C-C and C-H bonds (Dulev & Wi lliams! fl98ll) . 
Furthermore, PAHs and their cations were also supposed to ac- 
count for a subset of the "Diffuse Interstellar Bands" (DIBs) 



dLeger & d'Hendecourtl 


1 19851 Ivan der Zwet & Allamandolal 


1985 


Crawford et al. 




19851). more than ~ 300 absorp- 



tion features ubiquitously observed in the near-UV, visi- 
ble and near-IR in the spectr a of reddened stars (see e. g. 
Ehrenfreund & Charnle\l2000l and references therein). 

The RR nebula, besides being one of the brightest sources 
of AIBs in the sky, also has the almost unique prope rty of dis- 
playing a subset of emission ban ds in the visible dSchmidt et alJ 
1980; Van Winckel et al. 2002) which appear to be associated 
with a subset of DIBs dScarrott et al.lll992l). The RR and the 
R CrB star V854 Cen dKameswara-Rao & Lamberllll993l) are 
the only known cases in which DIBs have ever been possibly 
observed in emission. This makes the tentative identification of 
small neutral PAHs in this same nebula a potentially very ex- 
citing discovery, with far-reaching consequences also on the 
long-standing mystery of the identification of DIBs. 
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However, while the analysis in IViih etalJ fcOOfl narrows 
the range of potential carriers of the observed BL to a small 
number of molecules, it cannot provide a definitive identifi- 
cation, which requires an independent test. Such a test can 
be obtained from the IR and phosphorescence emission spec- 
trum of the same mole c ules. The low-fre quency vibrational 
modes (Langhoff 1996; Zhang et al. 1996) and phosphores- 
cence bands JsalmasCasffllQetaL 120041) provide an unique 
signature and ough t to be prod uced together with electronic flu- 
orescence bands (Brechignac 2005). 
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We carried out simulations of the photophysics of th e can- 
didate molecules using a Mon te-Carlo code llMulasiri 998: 
IJoblin et aljEool iMalloci et al.ll2003t iMulas et alJl2003l) . to- 
gether with quantum-chem ical calculat i ons for the r e levant 
molecular paramet ers JLanghoffl Il996t iMartin et alJ 1 19961 
Mall oci et all2004 and available laboratory measurements for 
the photoabsorption spectra Jjoblin et alJll992t Ijoblinl 1 1 992h 
and fo r visible and IR fluorescence quantum yields (Brech ignad 
2005). This produces a quantitative prediction of the IR and 
phosphorescence emission spectra for each given molecule, 
which must be related to the integrated BL attributed to this 
same molecule. 



2. Modelling procedure 

In previous pape rs Jjoblin et all 120021: IMalloci et alJ l2003t 
lMiilasetal.1 12003) we demonstrated the use of Monte-Carlo 
models to simulate the detailed photophysics of a specific 
PAH embedded in a given radiation field (RF). We here ap- 
ply the same procedure to the specific case of the three small, 
neutral PAHs which were tenta t ively i dentified in the Red 
Rectangle nebula by et alJ d2005b - namely anthracene 
(C14H10), phenanthrene (C14H10) and pyrene (C16H10) to de- 
rive their expected complete IR emission spectra. 

For the modelling procedure we used quantum-chemical 
results for structural paramet ers and vi brational analy sis, both 
from the literature (Langhoff 1996; Martin et al. 1996) and ob- 
tained by ourselves. For anthracene and pyrene we used th e ex- 
perimental photo-a bsorption cross-sections from Ijoblin et"afl 
( 1992); Joblin ( 1992J), w hile f or phenan threne we used the the- 
oretical spectrum from Mallo ci i et. al J |2004). We performed 
all of the quantum-chemical calculations in the framework of 
the Density Functional Theory (DFT), using the Octopus and 
NWChem computationa l codes. Deta ils on them can be found 
elsewhere (IMalloci et all2004l2005l) . 

A crucial modelling parameter is the assumed knowledge 
of the relaxation paths of the modelled molecules upon elec- 
tronic excitation following a photon absorption. Almost all 
neutral PAHs, including in particular anthracene, phenanthrene 
and pyrene, upon excitation S„ <— So, are well-known to un- 
dergo very fast internal conversion (IC) to the Si electronic 
level (see e. g. lLeachlll995l and references therein). From 
there, three relaxation paths are available, their relative weights 
dependent on the vibrational energy available: 

1. fluorescence Si — » So with a permitted electronic transi- 
tion, the remaining energy being subsequently radiated by 
vibrational transitions in So; relaxation of small PA Hs via 
this path is t he proposed origin of BL in the RR (HViih et alJ 
12004 1200.4 : 

2. direct Si "~> Ti or indirect Si T„ ~* Ti intersystem 
crossing (ISC), a radiationless transition followed by the 
emission of a phosphorescence photon in a Ti — > So 
spin-forbidden, electronic-permitted transition; the re- 
maining energy is radiated in vibrational transitions either 
(almost always) from So after the phosphorescence transi- 
tion or (very seldom) from Ti before it; 



3. internal conversion Si So, a radiationless transition, af- 
ter which essentially all excitation energy is radiated by vi- 
brational transitions. 

According to experimental results (Brechignac 2005J, the rate 
of fluorescence transitions (1) is essentially independent of the 
excitation energy. The rate of ISC radiationless transitions (2) 
increases slightly with excitation energy for the three molecules 
considered. The relaxation path (3) by IC to the ground state is 
open only when some excess vibrational energy in S 1 is avail- 
able, such threshold depending on the specific molecule and 
varying from ~ 2 10 3 crrT 1 for anthracene to ~ 4 10 4 cm -1 
for pyrene. Above this threshold the rate of IC to So (3) ex- 
ponentially increases, becoming by and large the dominant re- 
laxation path. The energy-dependent quantum yields for these 
thr ee relaxation paths were measured in gas-phase experiments 
by [Brechignac] J2005I) for all the three molecules under study 
here. 

All of the three relaxation paths described above con- 
tribute to IR emission from the ground So electronic state 
and are therefore considered in our model. Even the slow- 
est electronic transition, i. e. the Ti — > So involved in phos- 
phorescence, occurs wi th decay constants of the o rder of, at 
most, T ph ~ 1.5 10~ 3 s dSalinas C astillo et aTll2004l) . while the 
time constant involved in the fastest vibrational transitions is 
always tir > 5 10~ 2 s. This ensures that vibrational relax- 
ation almost always occurs after any electronic transition in 
the relaxation path, i. e. from So. Along any of the three re- 
laxation paths described above, the only process which can 
to some extent compete with the IR emission is the absorp- 
tion of another UV-visible photon before the end of the vi- 
brational cascade. When this happens, the residual excitation 
energy is added to that of the newly absorbed photon, mak- 
ing available higher-energy vibrational modes. Emission in the 
weakest IR bands, mostly towards the low-energy end, essen- 
tially happens only when the stronger ones, mostly towards the 
high-energy end, are energetically unaccessible; if vibrational 
cascades are interrupted before their end, weak bands are con- 
sequently suppressed, in favour of the stronger ones, which are 
emitted much more quickly when the molecule excitation en- 
ergy is high enough. 

We considered the modelled molecules in two different re- 
gions of the nebula and, therefore, our simulations correspond 
to two different exciting RFs. The first one is a Kurucz spec- 
trum with effective temperature T = 8250 K, sur face grav- 
ity lo gg = 1.5 and total luminosity of 6050 L Q rtVijh et alJ 
2005), summed with a blackbody at T = 60000 K to rep- 
resent the He white dwarf companion wi th a total luminos- 
ity of 100 L (Men'shchikov et al. 2002), both truncated at 
the Lyman limit and diluted ac cording to the geometry of 
the source llMen'shchikov et all 120021) at given angular dis- 
tances along the polar axes of the nebula. Extinction along 
the polar axes of the nebula is low and particularly gray 
dMen'shchikov et al-lEoo! IViih et aljEooll) and we neglect it 
altogether. We will therefore assume the spectrum of the es- 
timated RF to be constant along the bipolar lobes, and just 
scale with a wavelength-independent factor with distance, in 
the wavelength range which pumps the BL. 
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Table 1. Computed branching ratios for the three main relax- 
ation channels following UV/visible photon absorption by an- 
thracene, phenanthrene and pyrene in the RF of the lobes and 
halo of the RR nebula (see text for details). 



Relaxation branching ratios 
fluorescence I. S. C. I. C. 

(St-* S ) (S,^T,) (S,^S ) 



lomsation 







Anthracene (C 14 H 10 ) 




lobes 


0.157 


« 10~ 3 0.827 


0.016 


halo 


0.203 


« 10~ 3 0.796 


0.001 






Phenanthrene (Ci4H 10 ) 




lobes 


0.031 


0.372 0.587 


0.010 


halo 


0.042 


0.562 0.396 


3 10~ 5 






Pyrene (Ci 6 Hi ) 




lobes 


0.060 


0.284 0.642 


0.014 


halo 


0.094 


0.408 0.497 


0.001 



The second RF we co nsidered is the "attenuated" one, as 
shown in IViih et al.N200.4 . corresponding to what is seen by a 
molecule in the halo of the RR nebula, out of the bipolar cones 
and out of the dense, optically thick torus-like dust shell sur- 
rounding the central binary star. This RF is essentially due to 
heavily obscured light leaking through the torus and to light 
initially directed along the bipolar lobes and subsequently scat- 
tered into the halo. Since t he halo itself is optically very thin 
Jjylen'shchikov et afll2002t IViih et al"1l2005l) . we again assume 
the spectrum of this RF to be constant and just scale with a 
wavelength-independent factor with distance, over the energy 
range which pumps the BL. 

Both RFs are limited on the high-energy side by the Lyman 
limit at 13.6 eV, since higher energy photons are completely 
abso rbed in a small Hjt regio n close to the central evolved 
star (Men'shc hikov et al.l2002l) . Photons of energies lower than 
the absorption edge of the first permitted electronic transition, 
which for small, neutral PAHs is always higher than ~2 eV, are 
irrelevant for our purposes. The branching ratios for the three 
main relaxation channels following UV/visible photon absorp- 
tion by anthracene, phenanthrene and pyrene in the two RFs 
considered are listed in Tabled For the sake of completeness, 
we also included the ionisation yields, estimated according to 
the formula given bv lLePage etaP lEoOl). using the ionisation 
potentials 7.439 eV for anthracene, 7.891 eV for phenanthrene 
pyre n 

WebBook (Liasl2005l) . 

Under these conditions, the number 



and 7.426 eV for pyre ne taken from the online NIST chemistry 

^bi(v') 



dv' 



of BL photons 



radiated per unit time by a given molecule in the frequency in- 
terval between V and V + dv' will be proportional to the re- 
spective fluorescence quantum yield — fl *~ , multiplied by 

dv' 

the rate of absorption of exciting photons, in turn given by the 
product of the flux of exciting photons F exc (v) times the ab- 



sorption cross-section cr(v) of the molecule, i. e. 



dv 



= j dvF exc (v)cr(v) 
-- F exc J dv<& exc (v)cr(v) 



dv' 
dQ^v,v') 
dv' 



(1) 



exc dv' ' 



where we defined 
F 



F exc (v) 

dvF exc (v), Oe X c(v) = — 



and 



dv' ~ J 



dv<$> exc (v)cr(v) 



dv' 



With the assumptions we made for the RFs, F exc will vary with 
the position in the nebula, while O exc (v), which contains the 
frequency dependence of F exc (v), will only vary between the 
biconic lobes and the halo but will be position-independent 



within each of them. Since 



photon absorption rate, 



daft(v') 
dv' 



dv' 



is independent of the 



does not depend on the posi- 



tion along the line of sight. 

Completely analogous equations can be written for the IR 
emission and phosphorescence, i. e. 



^ir(v') - 



dv' 



with 



c/v<l> exc (v)cr(v) 



dQiR(v,v') 
dv' 



(2) 



— c/o^irO') 

F e . 



dv' 



^ir(v') 



dv' 



dv<& exc (v)o-(v) 



dQm(v, V) 
dv' 



and 

dR ph (v') 
dv' 



-I 



dv<& (v)cr(v)- 



dQ ph (v,v') 



_ — da ph (v') 

— t exc " 



dv' 



with 



dcrph(v') 
dv' 



I 



dv<b (v)cr(v)- 



dv' 



dQ ph (v,V) 



(3) 



~~dV 



As explained above, 



dQmjvy) 
dv' 



does slightly depend on 

the photon absorption rate if the average time between pho- 
ton absorptions is smaller than the average time it takes for 
an excited molecule to complete its IR emission cascade. 
dQphiv, V) 

This is not the case for , since phosphorescence, 

dv' 

although slower than fluorescence, still occurs on a much 
shorter timescale than the average interval between two pho- 
ton absorptions in all cases considered here. The detailed 
three-dimensional distribution of the molecules supposedly 
emitting the BL along a given line of sight is unknown: the 
emitting molecules might equally well be concentrated close to 
the source, where their photon absorption rate would be higher, 
or in a large volume extending relatively far out, with corre- 
spondingly lower photon absorption rates. We therefore con- 
sidered a range of distances to the central source, and hence of 
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photon absorption rates, evaluating the corresponding variation 



of 



d Qm{v, V) 
dv' 



and of the calculated IR emission spectra. 



Under optically thin conditions the observed fluorescence 
photon flux per unit solid angle on the sky along a given direc- 
tion is given by: 



dF u (V) 
dQ 



An J' 



drn(r) 

1 tfiffl(v') 



4-n dv 
1 afoW) 



- f 



d R bl (r, V) 
dv' 

drn(r)F exc (r) 



(4) 



An dv' 
with the definition 



T, 



T 



drn(r)F exc (r). 



Equation © assumes — ^ — - to be constant along the line 

dv' 

of sight, which is appropriate for a line of sight traversing only 
the halo of the RR. For a line of sight going through both the 
halo and one of the biconical lobes, the integral over the line 
of sight can be split in a part in the lobe and a part in the halo, 



in each of which 
generalised to 

dF u (V) = J_ 
dQ An 



dv' 



is constant. Equation © can be thus 



/ j — lobe / /\ 

I do= a (V) 



Ylobe + 



i — halo / / \ 

dan (v ) 



Thalo 



(6) 



dv' l " dv' 

where Ti b e and Yhaio are still defined by Eq. l|5}, the only dif- 
ference being the domain of integration, which includes respec- 
tively only the part of the line of sight in the lobe or only the 
part in the halo. 

Equivalent equations can be written again for both the IR 
emission and phosphorescence: 



dF lR (v') 
dQ 



1 

4^ 



dF ph (V) _ J_ 
dQ ~ An 



dv' 

P h (V) 



T lobe + 



dv' 



r h3 



( j — lobe 

dcr. 



dv' 



If lobe + 



j — halo / / \ 
<^ph (v). 



halo 
J 
\ 



dv' 



"Tfhalo 



(7) 



(8) 



In Eq. @ above 



do^ e (v') 



and 



do^R°(v') 



dv' dv' 
effective values, averaged along the line of sight 



are supposed to be 



The quantities 



and 



can be estimated 



dv' dv' 
for anthracene, pyrene and phenanthrene using the assumed 

spectrum of the exciting r adiatio n fields and the quantum yields 

measured by Brechignac] J2005t) . If we consider a line of sight 

dominated either by the lobes or by the halo, we can use the 

simpler Eq. ®. Integrating both sides of Eq. Q over V we 



obtain 



dF u {v') 
dQ. 



An J V dv' 

/■ 



(9) 



T 
4^ 
T 
4^ 
T_ 

An 



<ivO exc (v)o-( 



v) j dv' 



dQn(v,v') 
dv' 



I 



dv<& exc (v)o-(v)Qn(v) 



where the right hand side is the integrated BL photon flux due 
to a given molecule, observed along a given line of sight and 

cFfl is defined by the equation above. 

Th e BL spectrum in the R R has been reported bv lViih et alJ 
d2004h and lViih etail J20Q5I) . If a fraction 77 of the total, inte- 
grated BL photon flux observed is due to a given molecule, we 

T 

can solve Eq. with respect to — , to yield 



4tt 



i = ± r 

An cr fl J 



dv 



dF bi (v) 
dQ ' 



This result can be plugged back in Eqs. @ and ( 
we obtain respectively 



( 5 ) dF 1R (v') _ do'iRiv') T] 



dQ 



dv' cr fl 



dv 



dF ph (v') <io = p h(v') 7] 



(vl_n_ r 
J ^fl J 



dv 



dF h \{v) 
dQ 

dF h \(y) 



(10) 

\, from which 
(11) 



(12) 



Integrating Eq. fl!2i and proceeding as in Eq. we obtain: 



r^, dF ph (v') _ r da ph (v') t? r 

J V dQ ~ J V dv' v n J 



o" P h 



0"fl J 



dv' cr fl 
dF hl (v) 



dv 



dF bi (v) 
dQ 



(13) 



dv 



dQ 



The ratio , for a given molecule in a given RE is simply 
0-fi 

the ratio of the branching ratios for phosphorescence and flu- 



_ ^^(V') 

orescence, as listed in Table [T] The quantities 



1 — halo / / \ 

da'iR (v) 
dv' 



dv' 



and 



are part of the results of our Monte-Carlo model. 

The integrated BL photon flux along the lines of sight at offsets 
of 2.5"S 2.6"E and 2.5"S 7.8" E from the centra l so urce of the 
RR ne bula were obtained from lViih et alJ J2005I) and lViih et alJ 
assuming an average energy of the BL photons of 
~3.1 eV. They are respectively ~ 9.2 10 9 photons s _1 cm 2 sr _1 
at 2.5"S 2.6"E and ~ 2.5 10 9 photons s" 1 cm" 2 sr 1 at 
2.5"S 7.8"E. 

The line of sight at the offset 2.5"S 7.8"E is completely in 
the halo of the RR nebula. Therefore in this case, with the ex- 
ception of r], all quantities on the right hand side of Eqs. dl It 
and (I12> are known or can be calculated. We can thus quantita- 
tively estimate the expected IR emission and phosphorescence 
spectrum of anthracene, phenanthrene and pyrene apart from a 
scaling factor 77. 

The line of sight at the offset 2.5"S 2.6"E traverses both the 
southern lobe and t he halo of the RR nebula. Given the geom- 
etry of the source dMen'shchikov et allEool . the part in the 
halo of this latter line of sight is very similar to the line of sight 
at the offset 2.5"S 7.8"E; we hence assume their contributions 
to both the BL (~ 2.5 10 9 photons s~' cirT 2 sr~') and the IR 
emission to be approximately the same. The remaining BL at 
2.5"S 2.6"E (~ 6.7 10 9 photons s _1 cm 2 sr 1 ), after subtract- 
ing such estimated contribution from the halo, is due to the 
portion of the line of sight through the lobe. Therefore for this 
part we can again use Eq. (HTY with the RF in the lobes. The 
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total estimated IR emission spectra at 2.5"S 2.6"E can there- 
fore be obtained from the sum of the estimated spectra from 
the portions of the line of sight respectively through the lobe 
and through the halo. 



3. Results 

3.1. IR emission 

Tables[2E]E]and|5]list the most intense mid- and far-IR emis- 
sion bands expected, along with their calculated fluxes, respec- 
tively for anthracene (using two different vibrational analyses), 
phenanthrene and pyrene, assuming each of them, in turn, to 
be the sole carrier (i. e. rj =s 1) of the observed BL at po- 
sitions 2.5"S 2.6"E and 2.5"S 7.8"E from the central source. 
These fluxes can be simply scaled by the appropriate factor if 
the given molecule is assumed to only contribute a fraction 
of the luminescence. Since the IR fluxes scaled in this way 
still do slightly depend on the assumed photon absorption rate 
(see the discussion in the previous section), we ran independent 
Monte-Carlo simulations assuming a range of different photon 
absorption rates. In the halo the photon absorption rates turn 
out to be always small enough to allow almost all vibrational 
relaxation cascades to finish. This means that the estimated IR 
emission from the halo is independent of the three-dimensional 
distribution of the molecules along the line of sight. The situ- 
ation is slightly different for the lobes. For the closest position 
to the central source along the line of sight at 2.5"S 2.6"E we 
estimated an average time between absorptions r a b s ~ 70 s for 
pyrene and t r \, s ~ 60 s for the other two molecules. We com- 
puted the estimated IR emission fluxes for the portion in the 
lobe of the line of sight at 2.5"S 2.6"E assuming the above ab- 
sorption rates and those given by an RF dilution by a factor of 
~ 6, corresponding to molecules distributed farther away from 
the nebular axis, near the edges of the lobe. It turns out that only 
the very weakest IR-active bands of each molecule are signifi- 
cantly affected, being somewhat suppressed (up to a factor ~ 1 .6 
in the worst case) when using the higher photon absorption 
rates. For all the other IR-active bands, the resulting variation 
in the estimated fluxes is smaller than the uncertainty of the 
calculated vibrational transition intensities we used. To evalu- 
ate the impact of the latter uncertainty, we performed two sets 
of simulations for anthracene, using the DFT vibratio nal analy- 
ses re spectively at the B3LYP/4-31G level of theory jLan ghoff 
1 19961) for one and at the B3LYP/cc-pvdz dMartin et alJll996h 
for the other. The latter is very nearly the best currently feasi- 
ble for these molecules, while the former is well known to pro- 
vide a very good comp romise between accuracy and compu- 
tation al efficiency JLanghofdll996t iBausch licher & Lan ghoff 
119971) . The overall agreement between the two calculations is 
very good: while the two calculations distribute intensities in 
a slightly different way among nearby vibrational modes, their 
total is almost the same (see e. g. the in-plane C-H stretches 
between 3.25 and 3.29 yum). The differences on single band in- 
tensities and positions are comparable to the relative inaccuracy 
of each o f the two upon comparison with available experimen- 
tal data (Hudgins & Sandfor dll998h . 
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Table 2. Absolute integrated IR emission fluxes expected for 
the most intense IR-active bands of anthracene, calculated us- 
ing the vibrational analysis performed at the B3LYP/4-31G 
level of theory. They were estimated for two offsets from the 
central source, with two different RF dilution factors for the 
position on the lobe (see text for details). 



Bandpos. (fim) Integrated flux/77 (Wsr ' cm 2 ) 

2.5"S 2.6"E 2.5"S 7.8"E 

73% lobe + 27% halo halo 
high abs. rate low abs. rate 



3.25 


1.1 10~ 8 


1.0 10~ 8 


2.0 10- 


9 


3.26 


6.9 10- 9 


6.7 10~ 9 


1.3 10- 


9 


3.28 


2.4 10-" 


2.2 10-" 


4.3 10- 


12 


3.28 


1.8 10~ 9 


1.8 10 -9 


3.5 10- 


10 


3.29 


1.0 10~ 9 


1.0 10- 9 


1.9 10- 


10 


6.16 


1.5 10- 9 


1.5 10~ 9 


3.0 10- 


10 


6.51 


4.2 10-'° 


4.2 10-'° 


8.9 10- 


1 1 


6.86 


8.3 10" 10 


8.5 10" 10 


1.8 10- 


10 


6.86 


3.9 10" 10 


3.9 10- 10 


8.2 10- 


1 1 


7.22 


3.0 10-" 


2.9 10-" 


6.0 10- 


12 


7.43 


5.4 10-'° 


5.3 10- 10 


1.1 10- 


10 


7.60 


8.7 10 -10 


8.5 10" 10 


1.8 10- 


10 


7.85 


1.1 10~ 9 


1.1 10~ 9 


2.4 10- 


10 


8.55 


1.5 10-'° 


1.6 10-'° 


3.3 10- 


1 1 


8.62 


4.0 10" 10 


4.0 10-'° 


8.5 10- 


1 1 


8.67 


1.0 10~ 9 


1.0 10- 9 


2.1 10- 


10 


9.95 


4.8 10- 10 


4.8 10- 10 


1.1 10- 


10 


10.41 


8.6 10" 10 


8.7 10-'° 


1.9 10- 


10 


11.00 


1.9 10-'° 


1.9 10-'° 


4.0 10- 


1 1 


11.32 


6.3 10- 9 


6.3 10- 9 


1.4 10- 


9 


12.56 


2.9 10~ 12 


3.4 10~ 12 


2.8 10- 


12 


13.71 


6.2 10~ 9 


6.3 10-' 


1.4 10- 


9 


15.33 


1.1 io-'° 


1.1 10- 10 


2.6 10- 


1 1 


16.34 


5.0 10" 10 


5.1 10- 10 


1.1 10- 


10 


21.26 


8.9 10" 10 


9.0 10-'° 


2.0 10- 


10 


26.35 


2.1 10-" 


2.9 10-" 


1.7 10- 


1 1 


43.68 


9.1 10-" 


1.3 10-'° 


3.3 10- 


1 1 


110.23 


5.0 10 _u 


7.5 10" 11 


4.0 10- 


11 



The resulting IR emission fluxes are shown respectively in 
Tables 13 and |51 Upon comparison, they show the accuracy of 
the vibrational analyses to dominate the uncertainty in the esti- 
mated IR emission fluxes. 

3.2. Phosphorescence 

Table[6]lists the integrated phosphorescence fluxes expected for 
phenanthrene and pyrene, again assuming each of them, in turn, 
to be the sole carrier (i. e. 77 ^ 1) of the observed BL at positions 
2.5"S 2.6"E and 2.5"S 7.8"Efrom the central source. As for the 
IR fluxes, phosphorescence can be simply scaled by the appro- 
priate factor if fluorescence by the given molecule is assumed 
to only contribute a fraction of the BL. Since the branching ra- 
tios for phosphorescence is extremely small for anthracene in 
all cases, it is omitted in this section. The average phospho- 
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Table 3. Absolute integrated IR emission fluxes expected for 
the most intense IR-active bands of anthracene, calculated us- 
ing the vibratio nal analysis p erformed at the B3LYP/cc-pvdz 
level of theory (iMartin et alJll996l) . They were estimated for 
two offsets from the central source, with two different RF dilu- 
tion factors for the position on the lobe (see text for details). 



Band pos. (/im) 


Integrated flux/?/ (Wsr 


'cm" 2 ) 




2.5"S 2.6"E 


2.5"S 7.8"E 




73% lobe + 27% halo 


halo 




high abs. rate 


low abs. rate 




3.26 


8.1 10" 9 


8.1 10" 9 


1.5 10" 9 


3.27 


9.0 10" 9 


9.1 10" 9 


1.7 10" 9 


3.29 


1.7 IO" 9 


1.7 10" 9 


3.2 10~ 10 


3.29 


1.4 10" 9 


1.3 10" 9 


2.8 10" 10 


6.14 


1.0 10" 9 


1.1 10" 9 


2.2 10" 10 


6.50 


8.8 10" 10 


9.1 10" 10 


1.8 10" 10 


6.94 


1.7 10" 10 


1.8 10" 10 


3.7 10"" 


6.98 


1.7 10" 10 


1.8 10"'° 


3.7 10"" 


7.13 


3.4 10" 10 


3.4 10" 10 


7.2 10" 11 


7.34 


5.0 10" 10 


5.2 10" 10 


1.1 10 1U 


7.71 


9.9 10" 10 


1.0 io- 9 


2.1 10" 10 


7.94 


9.9 10- 10 


1.0 10" 9 


2.1 10" 10 


8.61 


3.1 10- 10 


3.3 IO" 10 


6.7 10"" 


8.73 


1.1 10-" 


1.1 io- 9 


2.3 IO" 10 


8.93 


3.1 10" 10 


3.1 IO -10 


6.8 10"" 


10.01 


8.6 10" 10 


9.0 IO" 10 


1.8 10" 10 


10.50 


5.6 10- 10 


5.8 10-'° 


1.2 10" 10 


11.15 


2.6 lO" 10 


2.8 IO" 10 


5.9 10"" 


11.35 


4.9 10-" 


5.2 10" 9 


1.1 IO" 9 


13.76 


6.1 10-" 


6.4 IO" 9 


1.3 IO" 9 


15.43 


1.1 io- 10 


1.3 IO" 10 


2.6 10"" 


16.56 


7.4 10" 10 


7.8 IO" 10 


1.6 IO" 10 


21.28 


1.1 IO" 9 


1.2 10" 9 


2.6 IO" 10 


43.29 


6.9 10-" 


1.1 io- 10 


3.0 10"" 


111.11 


4.9 10-" 


7.4 10-" 


3.9 10"" 



rescence photon energy can be estimated from the phosphores- 
cence spectra published by Salinas Castill oet al.l J200 4) to be 
about 2.4 eV for phenanthrene and 2.0 eV for pyrene respec- 
tively. Such spectra were taken under experimental conditions 
specifically chosen to enhance the phosphorescence yield, in 
solution, meaning that we can only use them as a guide for the 
spectral profile of the expected phosphorescence of the same 
molecules in the RR. Indeed, gas-phase phosphorescence spec- 
tra, in a collision-free environment, can be expected to arise 
from highly excited vibrational states and therefore have a dif- 
ferent width and vibronic structure (Brechianac 2005). We em- 
phas ize that we rel ied only on the gas-phase measurements by 
iBrechigriacl d2005l) for the branching ratios, using Eq. Jl 3i . 

4. Comparison with available ISO data 

The RR nebula has been observed several times with both 
the Short Wavelength Spectrometer (SWS) and the Long 
Wavelength Spectrometer (LWS) of ISO. The resulting spec- 



Table 4. Absolute integrated IR emission fluxes expected for 
the most intense IR-active bands of phenanthrene, calculated 
using the vibrational analysis performed at the B3LYP/4-31G 
level of theory. They were estimated for two offsets from the 
central source, with two different RF dilution factors for the 
position on the lobe (see text for details). 



Band pos. (fim) Integrated flux/77 (Wsr ' cm 2 ) 

2.5"S 2.6"E 2.5"S 7.8"E 

73% lobe + 27% halo halo 
high abs. rate low abs. rate 



3.23 


1.3 10" 


-8 


1.3 10" 


-8 


1.9 10 


-9 


3.24 


1.5 10" 


-8 


1.5 10" 


-8 


2.3 10" 


-9 


3.25 


1.2 10 


-9 


1.2 10" 


-9 


1.6 10" 


10 


3.26 


2.1 10" 


-8 


2.0 10" 


-8 


3.1 10 


-9 


3.26 


1.9 10 


-8 


1.8 10" 


-8 


2.8 10" 


-9 


3.27 


8.1 10- 


-9 


7.9 10" 


-9 


1.2 10" 


-9 


3.27 


2.8 10" 


1 1 


3.0 10" 


1 1 


3.0 10" 


12 


3.28 


1.1 10- 


•9 


1.1 10" 


-9 


1.7 10" 


10 


3.28 


2.1 10" 


-9 


2.1 10" 


-9 


3.4 10" 


10 


3.29 


2.2 10" 


10 


2.3 10" 


10 


3.8 10" 


1 1 


6.21 


4.1 10" 


10 


4.2 10" 


10 


6.7 10" 


1 1 


6.23 


2.9 10" 


10 


2.8 10" 


10 


4.7 10" 


1 1 


6.27 


2.3 10" 


-9 


2.2 10" 


-9 


3.7 10" 


10 


6.43 


6.3 10" 


12 


2.9 10" 


12 


7.7 10" 


13 


6.57 


9.1 10" 


10 


8.9 10" 


10 


1.4 10" 


10 


6.68 


3.6 10" 


-9 


3.6 10" 


-9 


6.1 10" 


10 


6.84 


7.2 10" 


-9 


7.1 10" 


-9 


1.2 10 


-9 


6.93 


1.8 10" 


-9 


1.8 10" 


-9 


3.1 10" 


10 


7.04 


3.3 10" 


10 


3.1 10" 


10 


5.4 10" 


1 1 


7.06 


4.7 10" 


11) 


4.7 10" 


10 


7.9 10" 


1 1 


7.45 


1.1 10" 


-9 


1.1 10" 


-9 


1.9 10" 


10 


7.49 


2.3 10" 


1 1 


2.4 10" 


1 1 


4.9 10" 


12 


7.70 


7.9 10" 


10 


7.9 10" 


10 


1.3 10" 


10 


7.76 


2.6 10" 


1 1 


3.0 10" 


1 1 


5.3 10" 


12 


8.00 


4.6 10" 


-9 


4.6 10" 


-9 


7.9 10" 


10 


8.16 


4.4 10" 


10 


4.2 10" 


10 


7.3 10" 


1 1 


8.31 


9.6 10" 


10 


9.8 10" 


10 


1.6 10" 


10 


8.46 


3.1 10" 


10 


3.3 10" 


10 


5.5 10" 


1 1 


8.52 


1.4 10" 


1 1 


1.4 10" 


1 1 


2.9 10" 


12 


8.60 


9.6 10" 


1 1 


1.1 10" 


10 


2.0 10" 


1 1 


8.71 


5.8 10" 


10 


6.1 10" 


10 


1.0 10" 


10 


9.15 


3.9 10" 


10 


4.0 10" 


10 


6.8 10" 


1 1 


9.63 


1.6 10" 


-9 


1.6 10" 


-9 


2.9 10" 


10 


9.66 


2.9 10" 


10 


2.9 10" 


10 


5.1 10" 


1 1 


10.01 


6.1 10- 


10 


6.2 10" 


10 


1.1 10" 


10 


10.12 


2.4 10" 


1 1 


3.4 10" 


1 1 


6.5 10" 


12 


10.53 


1.9 10" 


•9 


1.8 10" 


-9 


3.4 10" 


10 


11.48 


4.4 10" 


-9 


4.5 10" 


-9 


7.9 10" 


10 


11.49 


6.4 10" 


10 


6.2 10" 


10 


1.2 10" 


10 


12.05 


4.7 10" 


1 1 


6.5 10" 


1 1 


1.6 10" 


1 1 


12.24 


2.2 10" 


-8 


2.3 10" 


-8 


4.0 10" 


-9 


13.58 


2.6 10" 


-8 


2.6 10" 


-8 


4.8 10" 


-9 


13.95 


8.6 10" 


10 


9.1 10" 


10 


1.7 10" 


10 


13.98 


5.0 10" 


11) 


5.1 10" 


10 


9.7 10" 


1 1 


14.12 


5.2 10" 


1 1 


8.4 10" 


11 


2.5 10" 


1 1 


15.94 


1.4 10 


-9 


1.5 10" 


-9 


2.8 10" 


10 


18.19 


1.9 10" 


10 


2.5 10" 


10 


6.0 10" 


1 1 


19.99 


3.2 10" 


10 


3.8 10" 


10 


8.4 10" 


1 1 


20.07 


1.1 10" 


-9 


1.1 10" 


-9 


2.2 10" 


10 


22.75 


5.4 10" 


10 


6.1 10" 


10 


1.3 10" 


10 


23.25 


1.2 10 


-9 


1.3 10" 


-9 


2.6 10" 


10 


24.72 


2.1 10" 


10 


3.2 10" 


10 


8.3 10" 


1 1 


41.08 


2.1 10" 


10 


3.5 10" 


10 


1.3 10" 


10 


44.26 


5.5 10" 


10 


7.1 10" 


10 


1.7 10" 


10 
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Table 5. Absolute integrated IR emission fluxes expected for 
the most intense IR-active bands of pyrene, calculated using 
the vibrational analysis performed at the B3LYP/4-31G level 
of theory. They were estimated for two offsets from the central 
source, with two different RF dilution factors for the position 
on the lobe (see text for details). 



Band pos. (jim) Integrated flux/77 (Wsr 'cm 2 ) 

2.5"S 2.6"E 2.5"S 7.8"E 

73% lobe + 27% halo halo 
high abs. rate low abs. rate 



3.25 


2.0 10~ 8 


2.0 10~ 8 


3.0 10" 


9 


3.26 


2.0 10~ 8 


1.9 10~ 8 


2.9 10- 


9 


3.27 


5.4 10-" 


5.3 IO- 9 


8.0 10- 


10 


3.28 


8.3 10-' 2 


8.8 10-' 2 


1.3 10- 


1 2 


3.29 


7.9 10-'° 


7.7 10-'° 


1.1 10- 


10 


6.26 


1.6 IQ- 9 


1.6 IO- 9 


2.6 10- 


10 


6.31 


3.7 10- 9 


3.6 10-' 


6.0 10- 


id 


6.78 


1.1 IQ- 9 


1.1 10-' 


1.9 10- 


10 


6.92 


1.5 10-'° 


1.5 10-'° 


2.3 10- 


1 1 


7.01 


3.2 10-'° 


3.4 10 -10 


5.6 10- 


1 1 


7.01 


3.3 10" 9 


3.3 10~ 9 


5.6 10- 


10 


7.61 


1.8 10-" 


1.8 10- 9 


3.0 10- 


10 


7.98 


9.8 10-'° 


1.0 10~ 9 


1.7 10- 


10 


8.29 


3.0 10~' 2 


4.5 10~' 2 


1.1 10- 


12 


8.42 


2.8 10- 9 


2.8 10- 9 


4.7 10- 


10 


8.62 


4.3 10-'° 


4.2 10 -10 


7.4 10- 


1 1 


9.16 


1.2 10-" 


1.3 10~ 9 


2.2 10- 


10 


10.04 


1.4 10-'° 


1.5 10-'° 


2.6 10- 


1 1 


10.25 


6.5 10-'° 


6.5 10-'° 


1.1 10- 


11) 


10.47 


5.2 10~" 


5.5 10-" 


1.0 10- 


1 1 


11.79 


2.7 10~ 8 


2.7 10- 8 


4.8 10- 


9 


12.20 


6.8 10-'° 


7.0 10-'° 


1.3 10- 


10 


13.40 


2.0 10" 9 


2.0 10~ 9 


3.7 10- 


10 


14.06 


6.2 10-" 


6.3 10- 9 


1.1 10- 


9 


14.43 


3.9 10-" 


5.5 10-" 


1.3 10- 




18.20 


4.3 10-'° 


4.4 10-'° 


8.2 10- 




20.00 


4.1 10-'° 


4.3 10-'° 


8.1 10" 




20.38 


2.7 10-'° 


2.9 10-'° 


5.7 10- 




28.32 


2.2 10-'° 


2.7 10-'° 


5.5 10- 




47.80 


4.6 10-'° 


5.3 10-'° 


1.0 10- 


10 


101.60 


9.2 10"" 


1.3 10-'° 


8.0 10- 


1 1 



Table 6. Absolute integrated phosphorescence fluxes expected 
for phenanthrene and pyrene. They were estimated for two off- 
sets from the central source (see text for details). 



Molecule 


Integrated flux/^ (erg sr 


1 cm 2 ) 




2.5"S 2.6"E 2. 


5"S 7.8"E 




73% lobe + 27% halo 


halo 


Phenanthrene 


0.438 


0.130 


Pyrene 


0.136 


0.035 



tra are available through the online ISO database. All of these 
spectra were taken with an entrance slit which includes by and 
large the whole RR nebula, e. g. 33" x 20" for the SWS and a 
circular aperture of 84" diameter for LWS. Therefore, to com- 
pare our estimated fluxes with ISO observations we would need 
to calculate the predicted IR emission spectrum on a grid of 
points adequately sampling the aperture on the sky of the ap- 
propriate ISO instrument. This, in turn, would require BL mea- 
surements on such a grid, which are not available to date. To 
get a (rough) estimate of the BL distributi on, we inter polated 
the available BL measurements from lViih et alJ J2005I) with a 
thin plate spline (TPS). A TPS is a smooth function and passes 
through the original points f(Xi,yi) = Zi- It also has the prop- 
erty to become almost linear in the independent variables when 
away from the points used to define it. It resembles the more 
commonly used cubic (or more generally polynomial) splines, 
with the additional qualities of being smooth (and not only con- 
tinuous up to a given derivative order) and naturally multidi- 
mensional, hence well suited to provide a representation of a 
2D quantity as the BL surface brightness on the sky. A TPS is 
defined as 

1 " 

f(x,y) = fl + aix + a 2 y+-Yj h A log(?f), (14) 
with the constraints 

n n n 

Yj>i=Y* htXi= lL bm=0 ' (15) 

;'=1 i=l i 

where rf = (x - Xj) 2 + (y- yd 2 . We used the implementation of 
the TPS provide d by the In teractive Data Language (IDL) en- 
vironment. Since Viih et al. (2005) provided BL measurements 
on two long slits both south of the central source, in the interpo- 
lation we assumed the BL emission to be symmetric for inver- 
sion. In order to obtain a sensible asymptotic behaviour of the 
extrapolated BL, we built the TPS on the logarithm of the BL 
surface brightness, instead of the BL itself. Hence, in Eqs. dl4t 
and J15i x and y are the angular offsets from the central source 
of the RR, and z is the logarithm of the corresponding inte- 
grated BL. Figure ^ shows a surface plot of the resulting TPS 
which, despite the small number of data points from which it 
was obtained, clearly shows the geometry of the RR nebula. We 
were therefore able to integrate the TPS on the ISO apertures, 
to yield respectively ~ 1.40 10" 17 W cirr 2 for the 33" x 20" 
rectangular aperture of SWS and ~ 1.45 lO -17 W cm 2 for the 
circular aperture of diameter 84" of LWS. We remark that these 
integrated fluxes ought to be considered with caution, given the 
very sparse spatial sampling of the BL. Nonetheless, we used 
them to scale the calculated fluxes listed in tables |3 E] and [5] 
which were computed for specific positions in the RR nebula, 
to these ISO apertures, for comparison with available ISO data. 
In particular, to estimate the fluxes expected to be observed by 
SWS using the 33" x 20" aperture from the fluxes calculated 
for the RR lobe 2.5"S 2.6"E from the central source, we mul- 
tiplied the latter by - 1.40 lO -17 W cirr 2 and divided the re- 
sul t by ~ 4.58 10~ 9 W cm -2 sr (the BL intensity measured 
bv lViihetalj f2005) at this position in the RR). To ease the 
comparison, Figs.l2lto[T6lshow the resulting scaled calculated 
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1 1 1 1 1 1 1 

1 1 1 II II 




- 


- ,V J"' 









Wavelength (y^m) 



Fig. 3. Same as Fig.[2]in the wavelength range 5-15 //m. 



E— W off&et (mm center (afesac) 



Fig. 1. Spatial distribution of the surface brightn ess of the BL, 
interpolated from the values given bv lViih e t al. (2005) using a 
TPS (see text for details). East is left, North is up, contour lines 
are in 5% intervals. 




Wavelength (/im) 



Fig. 2. Comparison between the estimated IR emission spec- 
trum of anthracene (C14H10) and an ISO spectrum of the RR 
in the wavelength range 2.5^.0 /urn. Calculated spectra, under 
different assumptions (see text for details), are drawn in dashed 
(calculated from the fluxes in second column of Table |2j, 
dash-dotted (third column) 

and continuous (fourth column) lines, the continuum-subtracted ISO 
spectrum is shown as a dotted line. The central positions of expected 
anthracene bands are marked by ticks. 



spectra, along with the relevant spectra taken from the ISO data 
archive. 

For the SWS we chose the spectrum with the TDT number 
70201801, a low-resolution full-grating scan. For the LWS we 
chose the spectrum with the TDT number 70901203, a medium 
resolution spectrum. In both cases, a cubic spline was fitted to 



nirr" — ' 

~ ii5l 



Wavelength (/im) 



Fig. 4. Same as Fig.|5]in the wavelength range 13-30 yum. 



the continuum below the emission bands, and subtracted be- 
fore comparison. No further processing was performed on the 
ISO spectra, which therefore are by and large the automati- 
cally reduced spectra as obtained from the ISO data archive 
tool, just continuum-subtracted. In each figure we show the 
calculated IR emission bands from a given PAH, respectively 
scaled from the "high absorption rate" and "low absorption 
rate" 2.5"S 2.6"E columns (dashed and dash-dotted line), and 
scaled from the 2.5"S 7.8"E column (continuous line). The 
continuum-subtracted ISO spectrum is plotted as a dotted line. 
The positions of calculated bands are marked by ticks on each 
plot. Each expected band is plotted as a lorentzian curve (cor- 
respon ding to homogeneous band broadening, see IPech et all 
whose integrated flux is normalised to the expected 
scaled value. The widths of the lorentzians were arbitrarily cho- 
sen to match those of the observed bands in the same spectral 
range. In particular, we used a FWHM of 0.024 /mi for the 
in-plane C-H stretches around 3.3 /mi, of 0.2 /mi for bands 
between 6 and 15 /mi, of 0.3 /mi for bands between 15 and 
25 /im, of 0.5 /mi for bands between 25 and 50 /mi and of 0.9 
for bands longwards of 50 /mi. These values are just an edu- 
cated guess: many factors concur in determining the band shape 
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Wovelength (/^m) 




3.0 3.2 



3.6 3.8 



Wavelength (yLtm) 



Fig. 5. Same as Fig.[2]in the wavelength range 25^-5 yum. 



Fig. 7. Same as Fig.[2]for phenanthrene (C14H10) in the wave- 
length range 2.5^4.0 fim. 



a. 
'e 
s 




Wovelength (fJ-m) 



^ 5.0*1D" — 



Wavelength (/im) 



Fig. 6. Same as Fig.[5]in the wavelength range 100-120 jjm. 



Fig. 8. Same as Fig.0in the wavelength range 5-15 pm. 



and width for IR fluorescence bands of PAHs, such as lifetime 
broadening, anharmonicity, rotat ional structu re. For a discus- 
sion of these effects see e. g. lPech et al.l ( 120021) . Such a detailed 
treatment would only be feasible, in principle, for those spe- 
cific bands for which experimental measurements at different 
temperatures are available, i. e. not for the far-IR ones which 
appear to show the strongest diagnostic capability. Moreover 
it would be beyond the scope of the present work, since the 
precise band shape does not greatly affect its detectability. 

We just remark that, on general grounds, bands which 
are mostly emitted when molecules are more highly excited 
are expected to be slightly redshif ted due to anharmonicity 
dJoblin et al]ll995l IPech et al.ll2002l) . As an example, we per- 
formed a detailed band shape calculation for the in -plane C-H 
stretch of pyrene, using the experimental data from ljoblin et alJ 
\ 1995b . Figure shows the result, which can be compared 
with Fig.[21to see the impact of anharmonicity. Since the re- 
quired experimental data are not currently available for phenan- 
threne and anthracene, we assumed for them the same anhar- 
monicity parameters of pyrene. The resulting band profiles are 
shown in Figs.ll8landll9l 



Considering the uncertainties in the calculated fluxes, an- 
thracene appears to be compatible with all the ISO spectra, 
regardless of whether we use theoretical DFT ban d posi tions 
or the available gas phase positions measured by ICane et all 
(1997) for the comparison with observations. The far-IR band 
at ~ 1 10 /urn might easily be hidden in the noise in the auto- 
matically reduced spectra available in the online ISO database. 
If anthracene alone were responsible for the whole observed 
BL emission, about two thirds of the total observed flux in the 
3.3 //m band would be due to it as well. The slight position 
mismatch seen in Fig.|2]can be easily due to anharmonicity, as 
shown in Fig. The band in the latter figure appears to be 
slightly too broad, but since we used the anharmonicity param- 
eters of another molecule this is not definitive. 

Phenanthrene, on the other hand, can possibly account at 
most for a small fraction of the observed BL: its rather strong 
quartet out-of-plane C-H bend at 13.6 fim and its skele- 
tal bending mode at ~100 /im are not apparently detected in 
the ISO spectra. The C-H bend was measured in vapour at 
~770 K to be at ~13 . 7 um (unpublished data from the me a- 
surements of Ijoblinl i 19921) and Ijoblin et alJ i 1994 Il995h . 
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Fig. 9. Same as Fig.[7]in the wavelength range 13-30 jim. 
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Fig. 10. Same as Fig.0in the wavelength range 25-45 yum. 

Moreover, if phenanthrene alone were responsible for the BL, 
it ought to produce more than twice the total observed flux in 
the 3.3 //m band. Even allowing for the fact that DFT calcula- 
tions at the 4-3 1G level are kno wn to somewhat ove restimate 
the IR -activity of C-H modes iLanghofll 1 9961 iHudgins etail 
12001 1) , this is still unrealistic. Finally, to account for the posi- 
tion mismatch between the predicted and observed band posi- 
tions would require an anharmonicity band shift about an order 
of magnitude larger than the one obtained using the parameters 
of pyrene. Summing up, we conclude that fluorescence by neu- 
tral phenanthrene might contribute at most a small fraction of 
the observed BL. 

The case of pyrene is similar to that of phenanthrene, al- 
though slightly less clear-cut: it should also show a band at 
~100 //m which is undetected. This is one of the very few 
PAH s for which far -IR gas-phase measurements are available 
llZhangetal .1 ^996). and the experimental band origin of this 
fundamental vibration lies at ~105.3 yum, unfortunately in the 
most noisy part of the ISO LWS spectrum available online. The 
expected flux in this band would be sufficiently quenched by 
a high photon absorption rate to make it compatible with the 
observed spectrum, if pyrene were mainly concentrated in the 
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Fig. 11. Same as Fig.^in the wavelength range 90-1 10 fim. 
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Fig. 12. Same as Fig.|5]for pyrene (CigHio) in the wavelength 
range 2.5^.0 fim. 

lobes and relatively close to the central source. The expected 
out-of-plane C-H bend at -11.8 /im is a stronger constraint. 
This band was mea s ured in vapour at ~600 K to be at ~ 1 1 .9 /mi 
Jjoblin etalJ 11994 Il995h . and, as mentioned in the discus- 
sion for phenanthrene, DFT calculations at the 4-3 1G level 
are known to somewhat overestimate the IR-activity of C-H 
modes. Even allowing for all of these uncertainties, the ex- 
pected band at ~1 1.8 //m is at least a factor of ~5 too strong to 
be compatible with the observed spectrum. As to the in-plane 
C-H stretch, if pyrene alone were responsible for the BL, it 
would produce all of the total observed flux in the 3.3 fim band. 
However, the band position we obtain even with a detailed band 
profile modelling which takes into account anharmonicity ap- 
pears to be shifted bluewards of the observed peak. Summing 
up, we conclude that neutral pyrene might contribute at most 
-20% of the observed BL. 

5. Discussion and conclusions 

IVijh et alJ J2004 and lViih et all J20051) put rather stringent lim- 
its on the plausible neutral PAHs as candidates for the BL 
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Fig. 13. Same as Fig.[21in the wavelength range 5-15 fim. 



Fig. 15. Same as Fig.^]in the wavelength range 25^5 fim. 



E 
a. 





Wavelength (/^.m) 



Wavelength (/im) 



Fig. 14. Same as Fig.[^]in the wavelength range 13-30 /mi. Fig. 16. Same as Fig.^]in the wavelength range 90-1 10 //m. 



in the RR: they restricted the range of possible candidates to 
molecules including from 3 to 4 aromatic rings. Among them, 
the only ones showing fluorescence spectral profiles compati- 
ble with the observed BL seem to be anthracene, phenanthrene 
and pyrene. 

Upon a quick examination of the browsable online data 
from the ISO spectral database and comparison with our es- 
timated IR fluxes, some bands stand out as the most promis- 
ing diagnostics, due to their relatively high estimated intensi- 
ties compared with ISO sensitivity: the out of plane C-H bends 
and the skeletal bending modes between ~ 1 00 and ~ 1 1 /urn. In 
particular, even a cursory comparison between the estimated IR 
fluxes for phenanthrene and the available ISO SWS and LWS 
spectra show that its contribution to the production of the BL 
must be very small: indeed, fluorescence from phenanthrene 
would be so inefficient both in the lobes and in the halo of 
the RR that huge column densities would be needed to yield 
a significant contribution to the BL, and they would in turn 
also yield IR band intensities incompatible with observations. 
Phenanthrene and pyrene would produce prominent phospho- 
rescence bands too, which have not been observed in the RR 
thus far d"Witdl2005h . All these elements, taken together, mean 



that only anthracene, among the three molecules considered, 
appears compatible with the observations. 

Anthracene, if responsible for the BL, would also produce 
a quite substantial fraction of the observed flux in the in-plane 
C-H stretch bands. If it were to produce the observed BL, it 
would at the same time produce about half of the total observed 
flux at 3.3 fim. This is in good agreement with the observed spa- 
tial c orrelation between the latter band and the BL dViih et alJ 
2005). 

However, we should also consider the bands which are ab- 
sent from the calculated IR emission spectra: with the excep- 
tion of the in-plane C-H stretch bands at 3.3 fim, all of the 
other classical aromatic bands, which are strong in the ISO 
data, are almost negligible in the calculated spectra for the three 
neutral PAHs considered here. If anthracene were responsible 
for the BL, then it would almost completely account for the ob- 
served flux in the 3.3 fim band; it follows that in this case the 
remaining aromatic bands must be produced by different carri- 
ers, which in turn must not contribute much flux to the 3.3 fim 
band. In particular, any mixture of neutral PAHs producing 
the observed flux in the out of plane C-H bends at — 11.3 fim 
would also necessarily produce a substantial contribution to the 
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Fig. 17. Comparison between the estimated IR emission spec- 
trum of pyrene (C16H10) and an ISO spectrum of the RR in the 
wavelength range 2.5^4.0 fim, including detailed band-shape 
modelling. Calculated spectra, under different assumptions (see 
text for details), are drawn in dashed (calculated from the fluxes 
in second column of Table 0, dash-dotted (third column) and 
continuous (fourth column) lines, the continuum-subtracted 
ISO spectrum is shown as a dotted line. The central positions 
of fundamental pyrene bands are marked by ticks. 
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Fig. 18. Same as Fig.[n]for anthracene (C14H10). 

3.3 /mi band, unless only large species, with more than ~40 C 
atoms, are present. A possible solution may be to consider PAH 
cations, whose in-plane C-H stretch bands are well-known to 
be much less intense with respect t o the same bands in their 
paren t neutrals (see e. g. iLanghofill 199(1 lAllamandola et all 
Il999i) . This might account for the other strong bands observed 
in the RR, including the 11.3 fim band, which is not as reduced 

by ionisation. 

Th is is consistent with observations by Bregma n et alJ 
i 1993b . which show the 3.3 and 11.3 fim bands to have very 
different spatial distributions, and would imply a scenario 
in which anthracene dominates the population of small neu- 
tral PAHs, producing both the 3.3 fim band and the BL and 
PAH cations produce the other aromatic bands, possibly with 



Fig. 19. Same as Fig.ll7lfor phenanthrene (C14H10). 



some contribution by large neutral PAHs for the 11.3 fim 
band. Some detailed modelling would be useful to determine 
whether such a scenario can be realistic and compatible with 
the photo-chemical evolution of PAHs in the RR environment. 

The present data provide a powerful, quantitative diagnos- 
tic which may be used as a definitive cross-check for the hy- 
pothesis that small, neutral PAHs are responsible for the BL 
observed in the RR. Firm conclusions will require a thorough 
and careful review of available ISO spectroscopic data. In par- 
ticular, all available LWS spectra in the 90-120 fim wave- 
length range ought to be appropriately reduced and stacked, 
to maximise the S/N ratio. The expected IR fluxes ought 
to be well within reach with the sensitivity of ISO obser- 
vations, and will definitely be observable with forthcoming 
Herschel observations, which will go two orders of magni- 
tude deeper (e. g. with the PACS instrument in the spec- 
tral range 57-210 fim, see the official PACS web page at 
http : //pacs . ster . kuleuven . ac . be). 

The accuracy of the IR emission spectra calculated in this 
paper would be greatly improved by a better understanding of 
the spatial distribution of PAH emission in the RR and, con- 
sequently, a better knowledge of the RF illuminat ing them. 
Indee d, the best available model of the RR Jlvlen'shchikov et all 
2002§ does not explicitly consider PAHs. Both more observa- 
tional work, to better sample the spatial distribution of the BL 
as well as that of the IR emission, and more modelling are 
called for. In particular, observations of the BL with a much 
high er spatial resol ution and sampling of the RR than those 
by tYljh et a fl J200.4 would be needed: this would enable us 
to adequately assess the dilution effect due to the large ISO 
S WS/LWS apertures without recurring to the extrapolations we 
used in Sect.|4] 

We emphasize that experimental measurements of the pre- 
cise position and intensities of all IR-active bands of an- 
thracene, phenanthrene and pyrene in astrophysically rel- 
evant conditions would be very useful and reduce mod- 
elling errors: while the spectral range correspond ing to 
"class i cal" AIBs has been well e xplored (see e. g. Upblin 
Il992l ISzczepanski & Vabl Il993t Lloblin et . alJ Il994l 1199.4 



Mulas et al.: Absolute IR and phosphorescence fluxes from PAHs in the Red Rectangle 



13 



Moutou et all Il996t ICook et alJ Il996l iHudgins etaD I l994| 



Hud gins & Allamandola 


1995 


Hudgins & Sandford 


1998; 


Allamandola et al. 


1999; 


Kim & Savkallv 2002: Oomens et alJ 



2003), gas-phase measurements of far-IR bands of PAHs are 
sorely lacking to date. Last but not least, gas-phase measure- 
ments of the phosphorescence spectra of phenanthrene, pyrene 
and possibly other PAHs, which are also lacking to date, would 
provide yet another independent constraint amenable to direct 
observational verification of their presence in space. 
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